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Abstract 

The nucleation process under different external conditions is con- 
sidered. It is shown that the duration of this process can be connected 
with the microscopic corrections to the free energy of the critical em- 
bryo. Connection between variations in the value of the critical em- 
bryo free energy and the duration of the nucleation stage is given for 
several types of external conditions. This connection is in some cases 
reciprocal to uncertainty relation in quantum theory. In Appendix the 
derivation of main features of the general theory on the base of restric- 
tions coming from the possibility of effective and stable observations 
is given. 

1 Introduction 

The first order phase transitions are usually studied at example of the tran- 
sition of the supersaturated vapor into a liquid state. This example allows 
to go away from the numerous parameters characterizing the state of the 
mother phase and the state of the new phase. But even in this case there is 
no perfect coincidence between theoretical predictions and results of exper- 
iments. Now it is clear that the stationary rate of nucleation is determined 
with a bad accuracy and there are serious physical reasons lying behind this 
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problem. So, one can speak about some uncertainty in determination of the 
stationary flow of embryos or of the stationary rate of nucleation. The bad 
accuracy can be caused by two main reasons 

• The absence of the real stationary conditions for nucleation. 

• The bad value of the free energy of the critical embryo. This value is 
included in the formula for the stationary rate of nucleation. 

It is clear that according to [IJ the transition of the embryos from the 
pre-critical zone to the post-critical can occur far from the position of the 
critical embryo and, thus, the formula for the stationary rate of nucleation 
has to be reconsidered. But this case is rather rare and here we shall consider 
situations when the transition from the pre-critical zone to the post-critical 
zone goes through the critical point, i.e. through the position of the critical 
embryo. 

The numerous investigations of the establishing of the stationary state 
in the near-critical region showed that the situation where the stationary 
state is not established are very rare also. To see such situations one has to 
cut off the power of metatstability practically immediately after the moment 
when this metastability was created. So, this opportunity is also out of 
consideration here. 

The paper is organized as following: In the next section the situation with 
the artificial cut-off of the supersaturation is analyzed. Then the situation of 
decay of metastable phase is considered. Here the form of relation connecting 
the variations of parameters of the process resembles the uncertainty relations 
in quantum mechanics. That's why in appendix the derivation of the basic 
characteristics of the calssical mechanics and quantum theory including the 
uncertaity relations is given. The last situation which is considered is the 
external conditions of graduate creation of metastability in the system. Here 
the form of relations connecting the same parameters of the process radically 
differs from those in the previous situation. 

2 Determination of the pure rate of nucle- 
ation 

The stationary rate of nucleation in the main order has rather transparent 
origin - the rate of nucleation is proportional to exponent of the free energy G 



2 



of the critical energy taken in thermal units. Here and later all values having 
the sense of energy are taken in thermal units kT, where k is the Bolztman's 
constant and T is the absolute temperature. 

Determination of the free energy G is a very complex procedure. The 
problem is that the critical embryo has a number of molecules i/c big enough 
to make useless all calculations based on dynamic laws of motion. On the 
other hand the number z/^ is not big enough to apply the laws of statistical 
mechanics. But since there is no alternative one has to use the approach of 
thermodynamic description. 

Unfortunately the situation is more dramatic because one has to calculate 
the exponent of the free energy. Although the relative error in determination 
of the free energy becomes small the exponent reflects the absolute errors 
and these errors are not small even with z/c — oo. 

Really, the extraction of a separate embryo from the whole system is 
some simplification. It works satisfactory because the intensity of exchange 
between the embryo and environment is much more weak than the intensity 
of relaxation in the embryo to the state of internal equilibrium 

^tch 

where tint is the time of internal relaxation and tch is the characteristic time 
of exchange with environment. But here appears a problem - it is impossible 
to determine concretely the type of conservation equations for the separation 
of the embryo. So, it is impossible to determine the ensemble in statistical 
mechanics. 

Prom the first point of view there is no problem because different ensem- 
bles in statistical mechanics give equivalent results. But this means only that 
the relative values of macroscopic variables are equivalent. More precisely 
the relative difference of values have the order of Inn/n where n is a number 
of particles in the system. Taking exponent one can see the difference in n 
times. 

Also one has to take into account the possibihty of fluctuations which 
gives for the thermodynamic potential the shift of the order n^/^/n i.e. 
So, the exponent will have the correction in exp(~ n^/^) times. This correc- 
tion is enormously big. Certainly, one can say that fluctuations are already 
taken into account but the trace of incompleteness of the theoretical deriva- 
tion still remains. 

The next source of inevitable difficulties is the limitations on the size of 
the embryo. There are two aspects of this problem 
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• Since the system has a finite dimension the continuous spectrum of the 
energy transforms into the discrete energy levels and all integrations 
have to be replaced by summations. It is vary hard to do because 
the Euler-Maclaurin's decomposition is not converging one, but only 
asymptotic and the formal transformation here is not possible. 

• There appear the simple geometric problems like the difference between 
position of the surface of tension and the equimolecular surface. The 
conception of the surface of the embryo which is necessary to write the 
term with the surface energy and the equation of the material balance. 
But positions of surfaces do not coincide. This leads to the additional 
term of the order of n}^^. 

Following the second item we shall write the formula for the free energy 

as 

G=-hp + av^l^ + cv^l^ 

where parameters a, 6, c have a simple physical meaning: h is the difference 
of chemical potentials in a mother and a new phase, a is the renormalized 
surface tension, c is connected with the difference between equimolecular 
surface and the surface of tension at the plane surface. Later there appear 
corrections due to the curvature of the surface and it is convenient to continue 
this decomposition writing it as 

oo 
i=-l 

The term with i = looks like const + Cq In z/ and it is connected with non- 
equivalence of ensembles. 

As the result of all these constructions one can state that there is a 
microscopic addition SG to the value of Go, there is also an addition Suc to 
the argument of maximum Uco- Here the values with a subscript are the 
values based on the capillary approach, i.e. on 

The experiments intended to get the rate of nucleation and the free energy 
of the critical embryo are ordinary constructed in a following manner: At the 
initial moment of time there is a metastable state with no embryos of the 
new phase. After some time tout the metastability in the system is artificially 
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diminished to forbid the formation of new embryos. Then the number of 
droplets will be 

Here no depletion of the vapor phase is taken into account. To neglect the 
depletion of the mother phase it is necessary to fulfill 

^cut ^ ^depl 

where tdepi is the time of depletion, which will be determined in the next 
section. Then here 

and 

SN = (exp(5G') - l)No 
For very small 5G one can linearize the exponent and get 

5N = 5GNo 



3 Decay of metastable state 

Now we shall consider the process of the mother phase depletion. To give 
quantitative estimates it is necessary to specify the rate of the droplets 
growth. For the supercritical embryos, i.e. for the droplets it is reasonable to 
adopt the free molecular regime of growth. Under this regime the question 
of profile of the mother phase around the droplet can be solved extremely 
simple - there is no such a profile and, thus, the mother phase depletion can 
be described in a very simple manner. 

Under the free molecular regime of growth the number of molecules inside 
the droplet grows as 

where ( is the supersaturation of the mother phase and r is some character- 
istic time which is approximately a constant value. 

Then the number of molecules inside the new phase will be 

Q= f Jit')vit')dt' 
Jo 

where 1/(1') is the number of molecules inside the droplet formed at t'. 
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The rate of nucleation is connected with the distribution function / over 
p = 1/1/3 by 

/ = Jr/C 

Then 

[z — xYdx 



Jo 



One can show that the depletion occurs in a very rapid avalanche manner. 
Before the essential depletion one can consider / as a constant /* at the 
beginning of the process and get 

Q = hz'/A 

The length of the spectrum is determined by the following condition 

^ r 

where 
Then 

and one can determine the time tdepi from the following condition 
where 

J* = /*C/t 

Having written J as Z exp{—G{i'c)) where Z is the Zel'dovich' factor one can 
come to 

.4 ,C,3 4C 



The last equation allows the analysis of variations of the time on depletion 
and the free energy of the critical embryo formation. 
Having inverted variations one can get 
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This is the final equation and one can see that the total number of droplets 
can be calculated as 

N — J^ftdepl 

In the main order 

N - exp(-3G/4) 

and 

5N ~ iVo(exp(-3(5G/4) - 1) 
The equation on SG, St can be linearized which gives 

This equation is reciprocal in its functional form to the uncertainty relation 
in quantum mechanics 

SESt — const 

which combines the uncertainty in energy E and in time t. That's why 
in Appendix the method based on uncertainty relation is presented. It is 
necessary to stress that the derivation in appendix has a special meaning 
and demonstrates some new features. 



4 Gradual creation of a supersaturation 

Ordinary the external conditions have a continuous slowly varying character. 
Then one can introduce the ideal supersaturation $, i.e. a supersaturation 
which would be in the system in the absence of formation of a new phase. 
Thus, the ideal supersaturation is fully governed by external conditions. At 
the variations of $ of a relative order of 

the behavior of $ can be linearized 

^t) ^ + —{t - t,) 
at 

Here * marks values at some characteristic moment. 
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The duration of the nucleation period is approximately 2tnuc where tnuc 
satisfies the evident relation 

f = t r-i 

''nuc I * S*-"- * 

So, the value of tnuc is absolutely independent on 5G and, thus, 

Stnuc 

The moment has to be chosen as the moment of the maximal intensity 
of the droplets formation, i.e. here as th moment of the maximal supersatu- 
ration. Then here 

dC 

-r = 
dt 

and 

dQ_M 
'dt ~~dt 

The variation 5G certainly exists, but does it take place the variation of 
G*? To see this variation one has to write the condition for the maximum of 
the supersaturation 



J,exp(r— = ^1 



One can get J* outside of the integral and see that the previous equation 
reduces to 

J* = slow function ~ const 

It means that J* is invariant. When we add to G some addition 5G nothing 
will be changed. Simply the moment will be attained earlier or later. Then 
it is possible to introduce a shift of and this shift will depend on 6G. 

As the result one can state that in this case there is no variations of 
parameters SG^, St. 



5 Conclusion 

Having analyzed three concrete situations one can see that only in the situa- 
tion of decay there is a variations of parameters and these parameters forms 
the relation reciprocal to uncertainty relation in quantum theory. 
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A The role of restrictions coming from the 
possibihty of stable and effective calcula- 
tions 

At first the aim of tliis review was to sliow tlie role of requirements coming 
from the possibihty of correct calculations in physics. Later the construc- 
tions based on the requirement to produce the stable calculations of the 
characteristics of the system gave some more general conclusions presented 
below. 

We start from the general point of view and instead of concrete phys- 
ical theory we consider the general qualitative causal theory based on the 
differential formalism. It is necessary to clarify the last sentence. When 
we mention "the differential formalism" it means that the theory uses the 
standard formalism of the differential calculation. The term "causal" shows 
that some events are considered as "reasons" of other events described as 
" sequences" . 

There is no other special assumptions to start our constructions, but 
later some rather evident notations will be made to give us the possibihty to 
present concrete results. 

A.l Causality and time 

It is necessary to introduce a variable (or a characteristic) to describe the 
property of causality. Really, our style of thinking is principally a causal 
style. But the qualitative formalism has to show what event is a reason and 
what event is a consequence. In everyday speech we use the terms "earlier" 
for the reasons and "later" for the consequences. For example, consider the 
number of cars and the number of crushes. It is clear and mathematical 
statistics can show that the number of crushes correlates with the number 
of cars. But what is the reason - the big number of cars or the big mimbcr 
of collisions? Mathematical statistics can not give us the answer. It attracts 
our attention only to the fact the the increase of cars is associated with the 
increase of coUisions. Certainly, we know that the increase of the cars is 
the reason of collisions (However, collisions are in some sense the source of 
beaten cars and leads to increase of the total amount of cars, but beaten 
cars are excluded from consideration). But how it can be proven? Only by 
the fact the the increase of the number of cars occurs slightly earlier than 
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the increase of the number of coUisions. In real macroscopic social systems 
one can not really observe this effect obviously and this produces additional 
difficulties. But our style of thinking is to search the reason and the reason 
is marked by the word " earlier" . 

It is quite possible to see the system where the number of collisions is 
the reason of the number of cars. Really consider the social or biological 
system where "collisions" are "sexual relations" and the role of cars is played 
by males and females. Then the opposite casual construction takes place. 
The number of collisions is the source of the number of males and females. 
And here one can say that collisions occurs earlier than the increase of the 
population. 

The straight result of the given example is the necessity to introduce the 
the characteristic responsible for the casuality. This characteristic is called 
as the "time". It will be marked by a letter t. The task of the theory which 
is going to be constructed is to determine the dependence of characteristic 
of the system x on t, i.e. x{t). If there are several characteristics X(i) of the 
system which are marked by the index i one can consider a vector x. At first 
we shall consider the case of one variable x. 

One can argue whether t is discrete or continuous. The arguments to 
consider discrete time t are connected with a quantum Zenon effect. To 
use advantages of the differential formalism we consider here t as continuous 
variable, at least at some first steps of our considerations. The interval for 
time will be [a, b]. Sometimes we shall take it as an interval [—1,1]. 

Here we have to state that the reversibility ordinary announced in the 
classical mechanics has absolutely another rather local sense. To see the for- 
mal reversibility one has to change all velocities and the there is no concrete 
method how to do this. So, this reversibility is only imaginary property. In 
a real world there is no way to change the direction of time. Time is the 
characteristic responsible for the causal relations in our world. 

A. 2 Properties of proximity 

The function x{t) has to be established by the theory and then it has to be 
checked by some experimental measurements. The measurement of x at the 
moment tt will be marked Xi. Certainly, there is a characteristic error 5t 
of the choice of the time moment. At the accurate measurements this error 
becomes infinitely small. Then to have an infinitely small error of x it is 
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necessary that x{t) has to be a continuous function 

X{t) e C[a,b] 

Here we have to choose the measure of proximity. According to the central 
limit theorem of the probability theory the distribution of errors of stochastic 
variable (let it be y) under some rather wide spread conditions goes to the 
normal distribution N 

N ~ exp{-a^{y - y f) 

where y is the mean value of y and a is some constant. 

The last relation leads to the choice as the most appropriate metrics the 
ordinary metrics 



\x\ 
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This metrics corresponds to the scalar product 

n 
i=l 

The last scalar product corresponds to the Pyphagorean theorem for or- 
thogonal basis 

where 

is treated as a vector sum and ttj, are some constants. It is necessary 
to stress that x^Jj^ will be orthogonal to all X(k) with k i,j and the given 
property can be used again and again. 

The necessity to use this property for our construction is the following. 
We have to stress that we do not know the "true" characteristics of the 
system. We can miss some of them. There is possible to see the situation 
when instead of a pair coordinates we take one coordinate which is a linear 
combination of the initial ones. But the form of the normal distribution N 
has to be the same as it stated by the central limit theorem. It is possible 
only when we take the mentioned scalar product. 

Really, for the probability of two independent characteristics and X(j^ 
we have 

P = PiPj = expi-ajxf^) exp{-ajxlj-j) = exp(-X(y)) 
Here we count from their mean values. 
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A. 3 Dimensionality of a physical space 

On the base of measurements we have to reconstruct the function x{t). Con- 
sider the simplest case which is the case of one material point, i.e. a simplest 
system without any external ad internal parameters and characteristics. How 
many coordinates is necessary to introduce in order to describe this system? 
The evident answer is that the simplest case is one coordinate. But this 
answer has one disadvantage which will be considered below. 

Suppose that x(t) is some signal which is governed by stochastic process of 
random motion. The results of Poia [2| show that for the stochastic walking 
the return back to origin with the probability 1 will be infinitely many times 
when the dimension of space disd = loTd = 2. When d = 3 or greater 
then the probability of the infinite number of returns is 0. 

The illustration of these results can be easily seen if we mention that the 
diffusion equation corresponding to this process has the Green function with 
essential part 

G'~exp(-^4)/4D(,)t) 

i 

where is corresponding diffusion coefficient. The rest in G is the nor- 
malizing factor depending only on t. 

We see that the functional form of G does not depend on the number of 
the spatial variables. Here lies one of the possible reasons why we see the 
diffusion process clearly. This functional form coincides with the functional 
form of the normal distribution. Again this form is the exponential of the 
square form of the variables. This allows to speak about the distributions of 
this form as the result of the random walking process and the fundamental 
functional form which will be used below. Again one can see the invariant 
character towards the choice of the variables or their ignorance. 

Now we have to describe the consequences of these results for the problem 
under consideration. When d = 1;2 the infinite number of returns allows to 
construct the infinite set of measurements at the moments of these returns 
and have all measurements as the zero values. So, the trace of the random 
walks disappears. This effect is unsatisfactory and we need to have d = 3 or 
greater to exclude this effect. 

Since there is no other characteristics of the material point there is no 
other candidates for the true dimension of the space and we have to admit 
that namely = 3 is the crucial dimension. 
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A. 4 Reconstruction of functional dependencies 

The number of measurements of trajectory x{t) is a big finite number going 
to infinity. On the base of these measurements one has to reconstruct the 
functional form x{t) (here it is sufficient to consider one variable). Since there 
are only two arithmetical operations (addition and multiplication) one can 
not go outside polynomials. Actually only polynomials can be constructed 
and calculated. All other functions which are ordinary used like sin, cos are 
no more than idealized infinite series of polynomials. 

We have to restrict the class of functional dependencies allowed for x(t). 
Really, for discontinuous functions one has to measure x{t) at every point 
t which is certainly impossible. So, the consideration of the class of con- 
tinuous functions is preferable not only from the physical point of view but 
also from the enormous expenditures of measuring. Fortunately according to 
the Weierstrass theorem every continuous function / at [a, b] can be approx- 
imated by polynomial P 

or 

maxa<x<b\f - P\< e^O 

The simplest form of approximation is interpolation. The property of 
interpolation means that having measured n times at moments ti,i = l..n 
the function x{t) we get Xi = x{ti) and the polynomial L has a property 
L(tj) = Xi- Certainly, it is possible to construct the polynomial of a power 
n — 1 in a unique form. 

But here one faces with the "no go theorem" which states that for every 
manner of the choice of interpolation nodes there exists a function which 
can not be interpolated [3]. Namely, the difference between function / and 
interpolation polynomial L has the estimate 

||/-L||c= O(lnn) 

where n is the number of nodes. 

There exist recipes of Feier and Valle-Poussin [Sj which allow to approxi- 
mate / in C, but these recipes have no property of interpolation: there exist 
nodes where P{ti) ^ Xi. So, we came to a strange situation with a trajec- 
tory which does not satisfy the results of measurements. Some analogy takes 
place in quantum mechanics. 
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Having started from the classical point of view we have to require that 
all measurements have to be satisfied precisely. 

It is not a accidental coincidence that in the method of Valle-Poussin the 
approximated function does not coincide with the measurements at the half 
of points. This number can be hardly decreased because the weight used in 
this method is the optimal choice [3]. 

We come to a strange conclusion that at least some measurements can 
not precisely define the investigated dependence. This can be explained by 
impossibility of the fully precise determination in experiments of all possible 
characteristics of the system. In quantum mechanics this effect is called as 
"uncertainty relations". 

A. 5 Restriction of possible trajectories 

The possible evident answer to solve the problems leading to a classical me- 
chanics appeared from the "no go" theorem is to consider instead of the 
space C[a,h\ the space of functions with restricted first derivative. For such 
functions one can see [4J that for every t 

\x{t)-L{t)\ < 0{\nn/n) 

where n is the number of nodes (the number of measurements). Here the 
nodes are the Chebyshev's ones. The interval is [—1, 1]. So, here the inter- 
polation procedure approximates the real trajectory. Now the (infinite) set 
of measurements can give us the form of trajectory. 
There appeared two important consequences: 

• There appeared a new auxiliary characteristic of the system - the veloc- 
ity V or the momentum p. Now the description has to take into account 
this characteristic explicitly. 

• The velocity of trajectory is limited by some constant c. This corre- 
sponds to the requirement of the special theory of relativity. So, one can 
assume that the special theory of relativity goes from this very simple 
restriction of the class of trajectories in order to have the convergence 
in procedure of interpolation. 

As the special result we come to a conclusion that the trajectory and the 
first derivative of trajectory are the basic characteristics in description of the 
state of the system (or of the particle). 
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One has to stress that p or v can not be considered as the variables fully 
equivalent to coordinate x. In classical mechanics there exists a picture of 
Hamilton where x and p are formally considered as a pair of coordinates. 
But one has to remember that initially x and p have different senses and x is 
the main variable, while p is additional one. In consideration presented here 
it appeared as the characteristic only because the restriction of the class of 
trajectories. 

A. 6 The configuration space and the phase space 

At first it is necessary to recall that there exists a simple style to present the 
state of a complex system - the configuration space and the phase space. 

The configuration space is a space where n is a number of coordinates 
of all particles in a system. The state of the system is a point in configuration 
space. The coordinates in configuration space are orthogonal and this is one 
of essential features in the future analysis. This orthogonality is an evident 
consequence of the simple fact that if different particles in the system are 
independent then the space is simply reduced to the sum of two configuration 
subspaces for particles (Certainly, the configuration space for a free material 
point is R^). We assume this fact at least for negligibly interacting parts of 
the complex system. 

When we add the momenta as auxiliary characteristics of the system we 
get the space i?^" of all coordinates and momenta of the system. Here also 
orthogonality takes place. One can say that this is simply the property of a 
linear orthogonal space or one can seek something behind this fact. 

The problem which appears for every theory is that one does not know 
the number of coordinates of the system. The simplest structure of the phys- 
ical construction implies that we have a system with the given number of 
interacting balls. They are referred as "particles". Certainly, one can not 
state that these particles are the simplest systems, they also have to be con- 
sidered as complexes with rather complicated structure. In the field models 
the number of degrees of freedom is principally unknown. So, our approach 
has to allow the generalization to unknown coordinates, it has to manifest 
some invariance for the squeezing and for the developing of description. In 
approach to construct the field theories by continual integration procedure it 
is also necessary that the kernel has to be invariant for the arbitrary choice 
of the number of coordinates [3]. 

The subspace responsible for internal degrees of freedom has to be sep- 
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arated in some sense from the subspace responsible for external degrees of 
freedom. The best way to ensure this property is to consider these subspaces 
as orthogonal ones. The property of orthogonality is associated with the 
Pythagorean theorem 

c = a + 

and with Euclidean postulate that at the given point one can put only one 
parallel line to the given hne. This postulate leads to precise recipe resulting 
in the unique geometric operation. Then this postulate can be considered as 
reflection in the everyday life of intention to construct quantitative theory 
describing the world. 

What functions satisfy the requirements put here? It is easy to see that 
only the square forms of characteristics (coordinates) allow the operation of 
squeezing. Really, for 

with some constants a(j) if we miss some the functional form remain the 
same and is instead of X(^j^,X(^k) we take a linear combination xi — bjX(^j) + 
bkX(^k) the form will be the invariant also. 

Certainly, the square form can be reduced to the sums of squares by a 
linear transformation. 

A. 7 Classical motion 

The knowledge of the function x{t) means that we know the functional form 
for x{t) which can be written as 

H{x{t),t) = const 

for some function H, where t means that all preceding times are involved 
in description. The function H is called an "energy" or having written as a 
formula it is called the " hamiltonian" . Giving more detailed description at 
the current moment t we can indicate the derivatives of trajectory 

H{x,dx/dt,d'^x/dt'^, ...,t) — const 

Since we have established that only the value x{t) and the value of the first 
derivative dx{t)/dt are included into description then we have 

H{x,dx/dt,t) — const 
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The explicit dependence on time has to be excluded since all behavior of 
the system has to be invariant in respect to t — to where to is the time of 
preparation of initial state. Then we get 

H{x,dx/dt) = const 

For the dependence on v = dx/dt since as it has been mentioned f is a 
formal auxiliary characteristic we have to take the formal square functional 
dependence as it is prescribed above 

H{x,v) = U{x) + 

where U is some function and m/2 plays the role of a(j) 
The requirement H = const or dH/ dt = leads to 



mdv /dt = — 




which the Newton's second law of motion. The trajectory of motion x{t) will 
be called the classical trajectory. 

The same transformations can be done for vectors x and v = dx/dt. 

The function — ^ is ordinary described as the sum of forces. Fortunately 
this function has a rather simple form. 



A. 8 Statistical mechanics 

We have marked that the complete set of characteristics of the system can 
not be presented even for the simple systems. The account of the missed 
degrees of freedom or missed parameters can be made in a manner like it is 
done in statistical mechanics. Here one can use the derivation presented in 
[6] . Instead of precise characteristics Xi) , V{i) one has to use the distribution 
function p which is the probability P to find characteristic q in the interval 
[g, q + dq] divided by dq 

p{q)dq = P{q G [q,q + dq]) 

When for H one can see 

H = Hi + H2 
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where Hi is the hamiltonian of the first group of particles and H2 is the 
hamihonian of the second group of particles then two subsystems are inde- 
pendent. It is shown directly from dynamic equation of motion. Then the 
distribution function p{H) has to be the product p{Hi)p{H2) 

p{Hi + H2)=p{Hi)p{H2) 

The last functional equation has an evident solution 

p ~ exp{i3H) 

which determines the distribution function with a parameter /3, which is 
proportional to inverse temperature. 

Certainly, in foundations of the statistical mechanics it is necessary to see 
the hypothesis of Boltzman which states that the averaging over trajectory is 
equal to the averaging over ensemble. Rigorously speaking one has to prove 
this fact based on dynamic equations but it is very difficult to do. 

Since the distribution has to be normalized 

J dTp{H) = 1 

where / dT is the integration over all possible states, one can determine the 
normalizing constant b in expression 

p = 6exp(— 

Then 

b=[J dTexp{-f3H)]-^ 

One has to show that the missed coordinates do not change the functional 
form of the distribution function. Really, for H = Hi + aq^ we have 

p{H) = p{Hi)exp{-(3ax^) 

and p{H) differs from p{Hi) only by the factor which coincides in the func- 
tional form with a normalizing factor. 

The same derivations can be done for the arbitrary choice xi = hjX(j) + 
bkX(^k) and the factors of the normal distribution appear. 

The coordinates f (j) can be separated due to the form of H. This leads 
to the Maxwell distribution over velocities f6j 

p(t;(i)) ~ exp(-mf(2.)/2) 

Summarizing this section we introduce the recipe to take into account the 
missed degrees of freedom by transition 
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A. 9 Simple solutions 

The linear dependence of H does not lead to any progress. Now we shall 
present the solutions initiated by the simplest square form of H, i.e. H = 
ax'^+bv^. From the last section to have the convergence we see that constants 
a and b have to has the same sigh as const in H = const. We have three 
possibilities 

• a = 

Then v = const and we have the straight uniform motion which is 
equivalent to the stationary state. 

• 6 = 

Then x — const and there is no motion. This case does not take place 

• a > 0, 6 > 0, const > 
The solution is 

X — A sin wt + B cos wt 
where w is the frequency of oscillations. 

The last case is the most common one. If instead of we choose Xj+i — Xj, 
where coordinates of particles we get the moving wave which satisfies 

the wave equation 

(9t^ dx"^ 

for a function / with a positive parameter k. 
A. 10 Uncertainty relations 

The base of the classical mechanics is in some sense contradictory. The 
value of velocity involved into the formalism of the classical mechanics is an 
ill defined value from the point of view of numerical methods. Really, the 
numerical definition of derivative 

, x(t + 5t) - x(t) 
V = lim ; 

(5t-»0 dt 
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has a problems in calculation ^(^+^'^)-^(^) fQj- small 5t. 

From a physical point of view of classical approach one can see that 
collisions with environment make the velocity of the Brownian particle a 
fluctuating value and the instant value of velocity can not be well determined 
while the coordinate of the Brownian particle is a stable characteristic. 

These features require to consider some characteristic values of fluctua- 
tions Sx and Sv of x and v. The concrete form of uncertainty relation can 
be established on the base of quantum mechanics which is derived in next 
sections. Here one can stop the narration of this section. 

These simple features show us that it is necessary to present some ap- 
proach which has to take into account the impossibility to have in one and 
the same moment the value of x and v. 

But also one can present some not so rigorous derivations leading to some 
interesting results, already the Liouville's theorem says that the volume in 
the phase space is conserved and, thus, SpSx which is the elementary volume 
is conserved. 

Analogous derivations following Hazen start from the Hamiltonian form 
of the law of motion 

dH _ , dH __ , 
dp ^ dx ^ 

where a sign ' marks the derivative on time and an evident Maxwell relation 
1^ = and one can see that ^ + ^ — 0. Since for every function 

opox oxap ox op 

f ^ — ^ one can find p'x" + x'p" = or ^{p'x') — 0. In terms of finite 
differences -^{SpSx) = states that SpSx remains some constant. 

Now one can consider relation SE5t. certainly, the system following dy- 
namic equations has 6E = 0. Then one can not use here Hamiltonian relation 
but simply consider E = E{x{t)). Then dE = ^^dt. Since |f = p' which 
can be regarded as a definition of momentum. Then dE — p'x'dt. Then 
5E5t — p'x'dtdt — p'dtx'dt — SpSx — h where h is Planck's constant divided 
by 2%. 

From the other point of view one can write Si^ST = 1, where is a 
frequency, T is a period. Certainly dT — dt and we see that SE — hSu. So, 
one can see that E — hu since the shift in potential, i.e. in E is possible. 
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A. 11 The amplitude of transition 

It is necessary to give a theory for an elementary system which does not have 
fixed X, V. Since x is regarded as a main characteristic and v as additional 
we imply that at ti there is xi and at t2 there is X2- It is necessary to get an 
amplitude of such transition. 

Since the system does not have fixed x, v one can not speak about the 
fixed point in the phase space but only about some probability to be in a 
fixed state in the phase space. We have already presented a transition from 
the fixed state to a distribution in a section devoted to statistical mechanics, 
but it is clear that this approach can not lead to a true result because in 
exp(— the main role is played by the states with a minimal H. The 
less is the energy the greater is the weight of the state. In construction of 
the distribution for a separate system one can not use H and has to replace 
it by the characteristic of the system which has minimum on the classical 
trajectory. Rigorously speaking this characteristic S has to manifest three 
properties 

• It attains minimum at the classical trajectory or somewhere near this 
trajectory 

• It has to be an additive function of two non interacting parts of the 
whole system 

S[1,2] = S[1] + S[2] 

• It has to be an additive function of time 

The last two properties are necessary to consider exp(a;S') as some ele- 
mentary probability or something connected with the probability. The first 
property ensures the correspondence between the classical theory and this 
approach. 

Fortunately it is easy to present such characteristic 

S{tQ,ti) = f\mv'^/2 - U{x))dt 

The Euler equation for this functional coincides with equation of motion. 
Then the amplitude K of transition can be presented as 

K{a, b) — ^ cxp(Q;S'(a, b)) 

All trajectories x{t) going from the state a to the state b 
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A. 12 Classical limit 



At first one has to see that the Euler equation for S is the classical equation 
of motion. The functional S can be written as 

rtb 

S = L{x, X , t)dt 

Jta 

where 

L = - U{x) 

Having considered a variation of trajectory x{t) x{t) + Sx{t) we get 

S(x + 5x) = S(x) + / [Sx'— + 5x—]dt 
J ox' ox 

To see this equation it is necessary to prove that x and x' are really indepen- 
dent variables but it i possible to do already in C°°. 
Then having integrated by parts and assuming that 

6x^^\ -Sx^^l -0 

which is evident since x is fixed at ta and tb one can see that 

d ,dL , dL 
dx 

and due to the arbitrary variation Sx one gets 

^dVdx'^ dx^ 

at the minimum (maximum) of S. Since L = mv'^/2 — U the last equation is 
reduced to 

dU 

mx = — — — 
ox 

which is the classical equation of motion. 

Now it is necessary to investigate the limit of this approach for macro- 
scopic systems to show that the trajectory goes to the classical hmit. 

The approximately additive character of potential energy U and kinetic 
energy ^iTnivf /2 is very important for future analysis. Then the energy H 
and action S are the additive functions also. 

For macroscopic objects the value of S is proportional to the number 
of particles N and, thus, the characteristic value of relative deviation Ax of 
trajectory from providing minimum is proportional to N~^/'^. It is very small 
and for macroscopic systems trajectory is the classical one. 



5S = -[5x[U^^' 
J ^dVdx" 
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A. 13 Functional integral 



The calculation of the amplitude 

K{a,b) = ^ exp{aS{a,b)) 

All trajectories x{t) going from the state a to the state b 

is rather difficult to fulfill and can be made by the formalism of the functional 
integration. Details can be found in [S]. For our purposes it is important 
that the integral 

K{a,b) = J exi){aS{a,b))D[x] 

where D[x] is the infinite number of differentials in every point of trajectory 
or in other notations 

/oo 
exp{aS{a, b))dx 
-oo 

is absolutely the same expression. 

The formalism of the functional integral is well defined only in some 
special cases. One can restrict the class of trajectories, for example, consider 
the broken lines. This way requires the special limitations on trajectories 
which are not known. But for the special subintegral functions one has no 
need to make such restrictions. Namely, for the subintegral functions of the 
gaussian form 

/ ~ expi^Square form of trajectory) 

one can define the functional integral. The cause is the Pyphagorean theorem 
which allows to squeeze the number of differentials. 

Fortunately, in the simplest cases of free dynamics (see the section " Sim- 
ple cases") the function S is the square function of coordinates and the 
functional integral can be calculated. 

One can see that the principle of squeezing coordinates is here the neces- 
sary condition to fulfill the calculations of the functional integral. Certainly 
here one has to observe a condition 

Re{a) < 

if S is restricted from below. 
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A. 14 Schrodinger equation 

At first one can see the principle of superposition 

K{b,a) = / K{b,c)K{c,a)dXc 

wliere tlie integration is taken over all positions of trajectory at tc- 

This property allows to introduce some basic amplitudes K{xi,ti;xo,to) 
giving them a name "the wave function" ip{xi,ti). Then instead of 

/■CO 
K{x2, t2] xi, ti)K{xi, ti] xo, to)dxi 
-oo 

we have 

ip{x2,t2) ^ / K{x2,t2;xi,ti)ip{xi,ti)dxi 

J —oo 

Consider S = J^^ L{x',x)dt for small intervals e = \ti — to I- Then 

^ 6 ' 2 ' 

where X is initial value of trajectory and y is the final value of trajectory. 
Then ^ 

ij{x, t + €) = j^^ — exp(eaL(^-— ^, ^^-^))^{y, t)dy 

and 1/Ais the normalizing factor. 

Taking into account the explicit form of L one can get 

/oo ( j^ — y)^ X ~\~ y 

^ -exp(a — '— )exp(ea[/(^— 

Having introduced rj — y — x one gets 

— exp(Q;— — ) ex^{eaU{x + 77/2, t))ip{x + rj, t)dr] 

-00 ^€ 

Now one can decompose if] in powers of e and get 
^{x,t)+e = J_^jexp{a^){l + aeU{x + ri/2,t))i;{x + ri,t)dri 

One can see that for 

2ea'm ~ 77^ 
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the characteristic cancellation takes place and, thus, rj has also some small- 
ness. So, it is necessary to decompose in powers of rj which gives 



and 

7/>(x,t)+e^^^^ = l^^^expia^){l+aeU{x,t))[i;{x,t)+r]^+^r]''^]dri 



In the zero order 



and then 



, , , f°° 1 . mrf , , , / , 
= -rGxp{a——)dr]tp{x,t) 
D A It 



A = { -) 

mot 



Having calculated integrals one gets 

dt dx'^ 

with l^l = |af/|, \B\ = \l/{2am)\ which is the Shrodinger equation. So, the 
dynamic equation of quantum mechanics is derived. 

This derivation reproduces the analysis presented in [S] but for an arbi- 
trary parameter a. Now the task is to determine a. 

The value of a can be determined by the fact that in squeezed dimensions 
the classical solution and the solution of Shrodinger equation must have one 
and the same form because there is absolutely no information what theory 
has to be applied. Classical and quantum approaches have to coincide in the 
parts where nothing is known about the system. 

The classical approach gives the solution in the case when nothing is 
known - this solution is described in the section " Simple solutions" . We know 
this solution - in classical approach this is the superposition of oscillations or 
waves. So, quantum approach has to lead to the same waves and oscillations, 
at least in their functional form. So, a has a purely imaginary magnitude 

a =~ i 
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A. 15 Sense of the wave function 



Having derived the dynamic equations one has to clarify the sense of the 
amphtude K or at least the wave function ■0. It is clear that these objects 
have to be connected with a probability. So, when we have quantum objects 
with the wave function ip in coordinate representation, i.e. as the function 
of coordinate x, we have the density n{x) as the probability Pdx to have a 
particle into interval [x,x + dx] as some function F oi ip 



One have to specify this function F and it can be done from some evident 
requirements coming from constructions of continuous (field) models. Cer- 
tainly, quantum approach has to allow constructions of continuous models 
at least to construct the field theory. Then it is necessary to fulfill the av- 
eraging over all wave functions ip. In statistical mechanics construction it is 
also necessary to average over all states, i.e. over all wave functions. (Here 
there is no necessity to consider restrictions to occupy one and the same en- 
ergetic level for different states providing different statistics, one can simply 
forget about them.) So, the formalism of the functional integration naturally 
appears here. The function to be averaged looks like exp(G) where G is pro- 
portional to H, S or some other similar function. For us it is only essential 
that there will be potential energy U which is included into these functions. 
Potential energy U has to be written via the density n as 



with some constant a. 

We do not know the real number of coordinates in D[iIj]. So, the required 
functional integration can be fulfilled only when the argument of exp has the 
gaussian form. Then it is necessary that 



The value of n has to be a real number which requires to have a real number 
in the rhs of the previous relation 



n{x) = F{i/j{x)) 




Then it is necessary to fulfill the functional integration 




n{x) ~ ip'^ 
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Alternative possibility is to take {Reipy which gives a hardly appropriate 
result of a quickly oscillating function. 

One can see that here again the Pyphagorean theorem plays the fun- 
damental role and the squeezing principle works to determine the sense of 
^. 

Having recalled that n is proportional to the probability we see that 
the wave function has a simple physical meaning: l^/'p is the differential 
probability to find a quantum object in an elementary interval near x. 

A. 16 Procedure of measuring 

Procedure of measuring is certainly the interaction between the quantum 
object under investigation and the macroscopic object giving the result of 
measuring. One has to realize that it is impossible to give the detailed the- 
ory of such interaction and certainly every type of measurements has its own 
particular features and, thus, the detailed theory describing the measure- 
ments. But one can come to the main features of the measuring procedure 
already from the general principles neglecting the concrete picture of inter- 
action between the quantum particle and the classical object. 

Let us speak about the quantum particle and the measuring device. The 
action of the device will be described by an operator A. The elementary 
interaction between the particle and the device will be presented by applica- 
tion of an operator A to a wave function ip. Since operator A is an arbitrary 
one, here no supposition is made. When we speak about the additive char- 
acter then it is reasonable to take ^4 as a hnear operator. The process of 
interaction between the particle and the device can not be controlled - one 
can not say how many times the particle interacts with the device. So, the 
final state has to be stable - the next application of A does not essentially 
violates the wave function. Then it is necessary that 

final ~ final 

Then if) final is an eigenfunction of A. 

The evident physical reason of such requirement is the observation of only 
resonances. It is known that the problem of eigenfunctions appears in clas- 
sical mechanics when the resonances of oscillating systems are investigated. 
Prom other side the resonances are the most clear observed peculiarities of the 
system and there is absolutely nothing strange that namely resonances are 
considered as the observable features while all other features are neglected. 
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It is known that a linear operator A has many eigenfunctions. What 
eigenfunction will be taken here? Now it is necessary to decide with what 
intensity (probability) the result of measuring is some ipi 

Aiijji = aii/ji 

where is an eigenvalue. 

We see that for the macroscopic flow of particles with wave functions ip 
(it will be marked as ^) the result will be ai functions ipi . Here o;, are some 
stable characteristics. 

After the procedure of measurement we have instead of ^ the sum J2i otii^i- 

i 

At first there is no requirement that there is an equality between ^ and 
'^iO^ijji. But if wc make a measTircmcnt by operator A of the characteristic 
a, then by operator B of characteristic h, etc, then it is necessary that the 
next measurement does not feel the previous one. For unique particle this is 
not true, but for macroscopic flow this has to be observed, certainly, at the 
imaginary level. We do not know who at when makes observations. But we 
have to adopt that somebody very curious makes this observation without 
any traces. Then we have to have to possibility to ignore this observation. 

The recipe of observation has to be one and the same for every observa- 
tion. It can not depend on the previous observation. 

To exclude the influence of the previous observation for the macroscopic 
flow we need a linear law of reconstruction \E' on a base of ■^j. Then ^ has 
to be reconstructed as a linear combination of -0^. Then 

* = ^ hii^i 

i 

and hi are coefficients in decomposition. Then since [■^p is the probability 
one can see that are the probabilities to get and the result of the 
measurement. Then = ai = hi. 

As the result we sec that the necessary features of observation can be 
estabUshed without detailed description of interaction. Certainly, to fulfill 
the requirements of real eigenvalues and completeness the operators A have 
to be self-adjoint ones in a corresponding space. 
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A. 17 Conclusions 



In this review all constructions are based on the conception of the measure- 
ments. At first it is introduced in classical mechanics. Then the quantum 

mechanics is regarded as some regularization of equations appeared in clas- 
sical mechanics. This way is used in other branches of science, certainly one 
can consider the equations of hydrodynamics as some very advanced way to 
make the equations of motion more stable. Here the same idea is used. 

All conclusions made above have to be checked many times before they 
can be regarded as reliable features of mechanics. But even now it is clear 
that attempts to build a primitive physical model for every phenomenon 
in nature and consider it as the absolute truth has an evident error in its 
foundation. 
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